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Abstract 

High-resolution, inelastic x-ray scattering measurements of the dynamic structure factor S{Q, uj) 
of liquid water have been performed for wave vectors Q between 4 and 30 nm~^ in distinctly 
different thermodynamic conditions (T= 263 - 420 K ; at, or close to, ambient pressure and at P 
= 2 kbar). In agreement with previous inelastic x-ray and neutron studies, the presence of two 
inelastic contributions (one dispersing with Q and the other almost non-dispersive) is confirmed. 
The study of their temperature- and Q-dependence provides strong support for a dynamics of liquid 
water controlled by the structural relaxation process. A viscoelastic analysis of the Q-dispersing 
mode, associated with the longitudinal dynamics, reveals that the sound velocity undergoes the 
complete transition from the adiabatic sound velocity (cq) (viscous limit) to the infinite frequency 
sound velocity (coo) (elastic limit). On decreasing Q, as the transition regime is approached from 
the elastic side, we observe a decrease of the intensity of the second, weakly dispersing feature, 
which completely disappears when the viscous regime is reached. These findings unambiguously 
identify the second excitation to be a signature of the transverse dynamics with a longitudinal 
symmetry component, which becomes visible in the S{Q,uj) as soon as the purely viscous regime 
is left. 

PACS numbers: 61.10.Eq, 61.20.-p, 78.70. Ck, 63.50. +x 
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I. INTRODUCTION 



Water occupies a prominent role in natural sciences, partly due to its relative abundance 
and central role for the existence of life on earth, but as well due to its many unusual prop- 
erties, which - despite intensive research efforts - still defy today a complete understanding 
It is believed that the peculiar physico-chemical behavior of water arises from its 
local structure which is characterized by an open three-dimensional hydrogen-bond network 
of water molecules with an almost perfect tetrahedral arrangement and a well-structured 
second coordination shell. The hydrogen bond is responsible for this highly ordered local 
structure, and it is therefore clear that the processes of H-bond breaking and formation play 
a central role in determining the dynamical and thermodynamical properties of liquid water. 
The lifetime of the H-bond is, at ambient conditions, in the picosecond range. Thus, the 
study of the water dynamics in the terahertz frequency range is of particular interest. Tradi- 
tionally, this is done by inelastic scattering techniques using neutrons (INS) or x-rays (IXS), 
or by molecular dynamics (MD) studies. More recently, the high-frequency dynamics has 
been studied experimentally by time-domain spectroscopy j^Q]. The key quantity in most 
of these studies is the so-called dynamical structure factor S{Q, u) {hQ and fkv = E denote 
the momentum and energy transfer, respectively), which is the time and space Fourier trans- 
form of the atomic density-density pair correlation function 5]. Starting from the pioneering 
work by Bosi et al. Ql and by Texeira et al. 0] , several neutron experiments ^, 0, and 
x-ray studies [ll, 12 .ll^l. 11^. 12 1 have been devoted to the study of the THz- frequency 



dynamics in liquid water. 

The most striking result of these INS, IXS and MD studies is the existence of a particu- 
larly large dispersion effect in the sound velocity as a function of frequency. More specifically, 
at T = 277 K and ambient pressure, the longitudinal sound velocity increases from its hy- 
drodynamic value Cq = 1500m/s to c^o ~ 3200m/s. This transition from the hydrodynamic, 
or zero-frequency, to the infinite-frequency sound regime is directly observed at Q = 2 nm~^ 
and energy E = 3 meV. This phenomenon bears strong resemblance with that observed in 
glass-forming liquids, where the sound velocity dispersion is due to the presence of a struc- 
tural (or a) relaxation process. If Tq, is the characteristic time of this process, the system 
has a solid-like elastic behavior for modes with frequency Q{Q) satisfying Q{Q)Ta » 1 and 
a viscous one for frequencies such that Q{Q)Ta « 1 (here fl{Q) = c{Q)Q is the excitation 
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frequency). A detailed viscoelastic analysis of the dynamic structure factor demonstrated 
that also in water the transition from low- to high-frequency sound velocity is actually driven 
by the structural relaxation process [3]. 

Temperature-dependent IXS studies, in which the density was kept approximately con- 
stant to 1 g/cm^ by changing the pressure, revealed that, as long as remains in the 
picosecond region (i. e. in the high temperature region), it follows an Arrhenius behavior 
with an activation energy comparable to the hydrogen bond energy. This result provided a 
link between the relaxation process and the hydrogen bond network, and offered an expla- 
nation for the microscopic origin of the structural relaxation process in water. On a time 
scale short with respect to the lifetime of these H-bonded local structures the collective dy- 
namics is very similar to that of the solid state, i.e. ice In the opposite limit, the local 
structures have sufficient time to relax, and therefore show a typical liquid behavior. In the 
intermediate region, for Q{Q)Ta ~ 1, the dynamics of the density fluctuations is strongly 
coupled with the making and breaking of the hydrogen bond network. 

A second important result, arising from the whole body of INS and IXS studies of liquid 
water, is the existence of a second excitation in the dynamic structure factor. This second 
excitation has a (almost) Q-independent energy, and becomes visible only at Q larger than 
^4 nm~^. This second feature, observed in the SiQ.uj) of ambient condition liquid water 
by INS 0, y, y, Q and IXsH Q, Q, has been attributed to transverse-like dynamics 
on the basis of a MD study [l^ where the comparison between longitudinal and transverse 
current spectra was performed. The existence of the second feature and its transverse 
origin can be explained within the same viscoelastic framework which has been successfully 
used to describe the longitudinal dynamics. Indeed, the transverse dynamics, known to be 
non-propagating in liquids, can be actually supported by the liquid structure as soon as 
the viscous regime is left and the response becomes solid- like. The "transverse" excitation 
acquires flnite intensity in the S{Q, u) - which is only sensitive to the longitudinal molecular 
motion - because of the lack of order in the structure, leading to a mixing of modes with 
longitudinal and transverse symmetry [3]. 

To summarize, the high frequency dynamics of liquid water seems to be controlled by a 
relaxation process -whose microscopic dynamics is related to the making and breaking of the 
hydrogen bond network- with a characteristic time r^. All the modes with frequencies such 
that Q{Q)Ta « 1 see a relaxing local structure (viscous regime), and behave liquid-hke: 
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low sound velocity and no transverse-like dynamics identifiable in the longitudinal current. 
On the contrary the modes with frequency satisfying Q{Q)Ta >> 1 see a frozen, solid-like, 
structure and behave as they were in the glassy phase: elastic response, high sound velocity 
and identifiable transverse-like dynamics. 

A different interpretation for the origin of the two excitations appearing in the S{Q,u}) 
of water has recently been proposed by Petrillo et al. |9] and by Sacchetti et al. who 
analyzed the INS and IXS data using a solid-like framework: the mode-mode interaction 
and the symmetry avoided crossing. In this context, the second dispersion-less feature is 
assigned to a local inter- molecular vibration, resembling one of the optic modes in ice jlflj ]. 
Moreover, at variance to Ref. 13,^^, where the longitudinal sound velocity dispersion is 
associated to the interaction of the longitudinal sound waves with a relaxation process, the 
transition from the hydrodynamic to the fast sound is attributed to the interaction between 
the sound waves and this optic-like mode. As the symmetry of the two modes is supposed 
to be the same, the two branches repel each other to avoid crossing, and the high frequency 
branch at large Q disperses with a slope larger than that of the low frequency branch at 
small Q. This model describes the INS and IXS spectra as well as the viscoelastic model; a 
definitive conclusion on the correct interpretation of the water dynamics can therefore not 
be drawn on the basis of a "best fit" procedure. Such a conclusion needs to rely on an ex- 
perimental basis, as that brought, for example, from the physical meaning of the parameters 
entering in the viscoelastic and in the solid-like models, and/or from their dependence on 
the thermodynamic state. 

In this paper we present new IXS data which significantly extend the thermodynamic re- 
gion investigated in previous IXS and INS experiments. Moreover, thanks to the increased 
performance of the IXS instrument, data with better statistical quality could be recorded, 
allowing an even more reliable extraction of the relevant parameters. The previous IXS 
study by Monaco et al. was focused on the low Q- region (mainly below 7 nm^^) where 
the second excitation does not give a significant contribution, and a consistent viscoelastic 
analysis of the whole spectrum could be performed. The present work covers a much larger 
Q-range (from 4 to 30 nm~^) and thermodynamic region by a choice of pressure and tem- 
perature such that the relaxation time changes by more than a decade. At variance with 
the previous works, this allows us to follow the evolution of the second weakly dispersing 
feature not only as a function of momentum transfer, but also all the way from the viscous 
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to the elastic limit. The extended thermodynamic and Q regions and the improved statistics 
allow us to provide overwhelming evidence in favor of a " viscous" - rather than solid-like - 
origin of the transition from "normal" to "fast" sound in liquid water and ii) the transverse 
origin of the low frequency weakly dispersing feature appearing in the S{Q,uj). The paper 
is organized as follows: In section II we present the experimental set-up and the information 
pertinent to the sample cell, the theoretical formalism utilized in the data analysis, and the 
experimental results. Section III provides the data analysis and the discussion, while Section 
IV draws the conclusions. 



II. EXPERIMENT 

A. Experimental setup 

The experiment was carried out at the Inelastic X-ray Scattering Beamline II (ID28) 
at the European Synchrotron Radiation Facility in Grenoble/France. The X-rays from an 
undulator source are monochromatized by a cryogenically cooled silicon (111) double crystal 
monochromator and a high-energy resolution backscattering monochromator, operating at 
the silicon (11 11 11) reflection order. The backscattered photons of energy 21.747 keV 
impinges on a gold-coated toroidal mirror, which provides a focal spot at the sample posi- 
tion of 270 (horizontal) and 80 (vertical) /im^ Full Width at Half Maximum (FWHM). The 
incident flux on the sample is 5x10^ photon/s. The scattered photons are energy-analyzed 
by a Rowland circle five-crystal spectrometer, operating at the same reflection order as the 
monochromator. The energy-analyzed photons are detected by a Peltier-cooled silicon diode 
detector which has an intrinsic energy resolution of 400 eV j|2(J]. The dark counts due to 
electronic and environmental noise amounts to about 0.003 counts/s. The momentum trans- 
fer Q = 2ki ■ sin{6s/2), where fcj is the incident photon wave vector and 6s is the scattering 
angle, is selected by rotating the spectrometer around a vertical axis passing through the 
scattering sample in the horizontal plane. Since there are five independent analyzer systems, 
spectra at five different momentum transfers can be recorded simultaneously. Their sepa- 
ration amounts to approximately 3 nm~^ for the Si (11 11 11) reflection. The energy scans 
are performed by varying the monochromator temperature while the analyzer temperature 
is kept fixed. Conversion from the temperature scale to the energy scale is accomplished by 
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the following relation: /S.E/E = aAT, where a = 2.58 • 10~^ is the linear thermal expansion 
coefficient of silicon at room temperature. The validity of this conversion has been checked 
by comparing the measured diamond dispersion curve for longitudinal acoustic and optical 
phonons with well established inelastic neutron scattering results. The overall experimen- 
tal resolution is experimentally determined by measuring the scattering from a disordered 
sample of plexiglass at a Q-transfer of 10 nm~^, corresponding to the first maximum in the 
static structure factor S{Q), and at T = 10 K in order to maximize the elastic contribution 
to the scattering. 

Distilled and deionized water was loaded into a specially designed stainless steel cell, 
made out of INCONEL-751, which can be pressurized by a hand pump up to 5 kbar and 
heated up to 700 K. The X-ray beam passes through two 1-mm-thick diamond single crystal 
windows with a 2.3 mm aperture. The distance between the two windows, i.e., the sample 
length along the X-ray beam was 10 mm. The cell geometry allowed us to cover the relevant 
momentum transfer regime up to 30 nm~^ with a sample length comparable to the X-ray 
photoabsorption length (yU ^0.1 mm~^ at E=21 keV). The high pressure cell was kept in 
vacuum in order to minimize both temperature gradients and air scattering. The pressure 
was monitored with a calibrated gauge whose precision is about 1%. The temperature was 
determined by a Cr-Al thermocouple in direct contact with the main body of the cell with 
an accuracy of ± 1 K. Table H] reports the thermodynamic conditions of the experiment as 
well as the respective shear viscosity (r/) and the adiabatic speed of sound (cq) values. 



B. Theoretical formalism and fitting procedure 

The theoretical framework used to interpret the experimental spectra is based on the 
generalized Langevin equation (memory function approach), more specifically on the gener- 
alized hydrodynamic description. In this framework the formal structure of the "classical" 
hydrodynamic regime is retained, but the thermodynamic derivatives and transport coeffi- 
cients are replacec 



[or frequency) 
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S functions which can vary both in space (or wave number) and time 
. In the last years, this approach has been successfully applied to 
the interpretation of the high frequency dynamics of different systems, including water, and 
spa„n,„g fr™, noble gases |24| to hqu.d ,„eta.s Q. as well as from supereoo.ed liquids Q 
to glasses j22l|. The basic equations are discussed in detail in Ref. here, we only stress 
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TABLE I: Thermodynamic conditions of temperature T and pressure P at which the IXS spectra 
were recorded. The uncertainty onf T is ±1 K and that of P is ±1 %. Columns 3 -5 report the 
three corresponding thermodynamic quantities used in the data analysis, namely the mass density 
p, the adiabatic sound velocity cq, and the shear viscosity r]. Their values were obtained from the 



water Equation of State 



2l| . The freezing point of water at P=2 kbar is 251 K. 



the key points and the details relevant to the present data analysis. 

The longitudinal dynamics is described by a Langevin equation for the density- density 
correlation function F{Q,t): 

^^^^ + u',iQ)F{Q, t) + 1^ dt'ruQit - t'f-^^^ = 0; (1) 

and, consequently, the expression for the S{Q,uj) (normalized to S{Q)), the time Fourier 
transform of F{Q, t), reads: 

S{Q,u;) ^ 1 ^oH(^) (0) 

where mQ(u;) and 771^(0;) are the real and imaginary part of the Fourier transform of the 
memory function mqlt). Its generalized form can be written as: 

mgit) = u^liQMQ) - l]e-^-(«)«'* + KLiQ,t), (3) 

where tu^ = KbTQ'^/ {MS{Q)) is the (squared) frequency of the sound waves excitation in the 
hydrodynamic regime, 7(Q) and Dt{Q) are the Q-dependent generalizations of 7 = Cp/C^ 
(the constant pressure to constant volume specific heat ratio) and Dt = K/{pCy), k is the 



thermal conductivity, S{Q) is the static structure factor, M is the molecule mass and Kb the 



Boltzmann constant. The second contribution to the memory 



directly related to the longitudinal kinematic viscosity 



'unction, KL{Q,t), a quantity 
, is here chosen as: 



I23|; 

Kr^{Q,t) = 2T^{Q)S{t) + A2(Q)e-*/-('3), (4) 

where 2Tf^{Q)S{t) represents a fast (microscopic or fi) decaying contribution to KiiQ^t) 
and the exponential term accounts for the structural (a) relaxation process. Tq, is the Q- 
dependent time which characterizes the long time tail of KL{Q,t), and A^(Q) is the struc- 
tural relaxation strength which is related to Coo(Q) and Co(Q) - the Q-dependent generaliza- 
tions of the usual infinite frequency and adiabatic sound speeds |28| via = (cqo — co)Q'^- 
In summary, the resulting expression for the memory function becomes: 



It is worth pointing out that the time dependence of the structural relaxation process contri- 
bution to the memory function is assumed here to be Debye-like (no stretching is present). 
This is in apparent contrast to recent findings P, where the existence of a non-negligible 
stretching {(3 ~ 0.6) in the correlation function of supercooled water was observed. However, 
as also pointed out in j^], the existence of stretching is hardly detectable in the tempera- 
ture region investigated here. As a matter of fact, we did not find any significant statistical 
improvement in introducing a new fitting parameter [(3). We therefore used a simplified 
version of the memory function and fixed (3 to one. 



The viscoelastic analysis of the water IXS spectra previously performed [15| was confined 
to Q- values not exceeding Q=7 nm~^, a region in which the second excitation discussed 
in the introduction is not visible. Consequently, the applied theoretical formalism yielded 
excellent agreement between the model function and the experimental data. In the present 
case, the inclusion of a second excitation is mandatory as it is evident, for example, by 
inspection of Figure Q which reports the IXS spectrum recorded at T=263 K, P=2 kbar 
and (5=10 nm~^. In this figure, panel (a) shows the best fit to the data with only one 
excitation (viscoelastic model) according to Eq. |21 Both the inelastic and the central part 
of the spectra are not well reproduced, and it is obvious that a second excitation has to be 
added. This has been done for the fit presented in panel (b) where the second excitation is 
accounted for by a damped-harmonic oscillator (DHO) function js^ added to the viscoelastic 
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FIG. 1: IXS spectrum at Q=10 nm~^, T=263 K and P=2 kbar. The experimental data and their error bars 
are shown together with the best fits utilizing the model function described in the text with and without 
the second inelastic contribution. Panel (a): one viscoelastic contribution. Panel (b): one viscoelastic 
contribution and one damped harmonic oscillator model. The inset shows the central peak region. 

S{Q,uj), resulting in a significantly better fit. It is important to stress that this is a purely 
empirical approach to introduce the mixing phenomenon of longitudinal and transverse 
dynamics. A correct theoretical description should provide a memory function, directly 
reproducing a double-excitation line shape, which naturally would involve a much larger 
set of variables. In particular, such a formalism should be extended to a regime of broken 
ergodicity, when the frequency range is much higher than the structural relaxation time 
{cuTa > 1). In this range the symmetry arguments that decouple the transverse (T) from the 
longitudinal (L) variables in the generalized Langevin equation must be abandoned and a 
coupling L — T mechanism must be taken into account. Work along these lines is currently 
in progress. Within the present context, the choice of a DHO lineshape to describe the 



9 



transverse dynamics is used empirically to provide information on the intensity and the 
frequency of the second excitation. In summary, the model function used to represent the 
spectra is composed of the following pieces: 

(1) A viscoelastic model function, proportional to Eq. |2l to account for the central peak 
and the longitudinal dynamics: 



We have neglected the contribution due to thermal relaxation, which amounts to set 
7((5) = 1, an approximation which turns out to be very good in the high Q region. 

(2) A DHO line-shape for the second excitation: 

Where Qt is the maximum of the "transverse-like" contribution to the the longitu- 
dinal current spectrum (the dynamic structure factor multiplied by tu^/Q^, i- e. of 
Lu"^ s^'^\Q , iu)) and Tt is the width of the inelastic transverse peaks. 

This model function for the dynamic structure factor is symmetric, and, in order to account 
for the quantized character of the energy transfers at the microscopic level, has to be weighted 
with a function that i) satisfies the detailed balance and ii) becomes unity in the classical 
limit {hjj /KbT << 1). Among the possible choices, the weighting factor usually utilized is: 

w{uj,T) = ^{n{uj,T) + 1) = - expi-hw/ksT)]-' 

Kb! Kb! 

Finally, the theoretical model has to be convoluted with the experimental resolution function 
R{uj), to give a fitting function of the form: 

I{Q,u;) = ® {^[n{u;,T) + l][s^'^\Q,u;) + s'^^\Q,u;)]} + B (5) 

Kb-I. 

Here B is an additional term which accounts for the electronic and environmental back- 
ground of the detectors. As a result, the data are fitted with 9 free parameters: i) the 
background B; the transverse ii) "intensity" At, iii) position Qt and iv) width Tt] v)the 
longitudinal intensity Al; the three parameters of the longitudinal memory function, namely 
the structural relaxation vi) time and vii) strength A and the area of the microscopic 
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relaxation F^; and, finally, ix) the static structure factor S{Q) entering in the term Uo- All 
these parameters are, in principle, T- and Q-dependent. 

As a matter of fact, F^, which represents the topological microscopic disorder contribution 
to the acoustic attenuation, is known to have a negligible temperature dependence. There- 
fore we have fixed its value to the one obtained at the lowest temperature. We checked that 
leaving this parameter free does not change significatively the values of the other parameters. 
Moreover, we observe that, for the highest temperature spectra, leaving the energy position 
and width of the transverse contribution as free parameters in the fit, results in a meaning- 
less overdamping of the DHO function to the extent that it cannot be distinguished from 
the purely relaxational elastic contribution already accounted for by the viscoelastic model 
(the intensity of this contribution vanishes at high T, thus implying that Qt and become 
irrelevant). In order to circumvent this, and guided by the observation of Sokolov et al. jssj ]. 
the position and width of the second excitation was unambiguously determined in an uncon- 
strained fit for the lowest temperature IXS spectra, and kept fixed at higher temperatures. 
Therefore, at all temperature but the lowest, there are six free fitting parameters. 

C. Experimental results 

IXS spectra were recorded at the different temperatures and pressures as indicated in 
TablelU spanning the momentum transfer, Q, region from 4 to 16 nm~^ with an approximate 
constant spacing of 3 nm~^. For the lowest (263 K) and highest (419 K) temperatures, the 
Q-range was extended to 30 nm~^. The spectra extend up to 60 meV on both sides of the 
elastic line, with an integration time of 100 s per spectrum. For each setting three spectra 
were collected, and subsequently summed in order to achieve the high statistical accuracy 
required for the data analysis. Typical total counts range between 1500 (at 4 nm~^), 2500 
(at 9.93 nm~^), and 3500 (at 12.91 nm^^). To account for the slow drift of the photon flux 
impinging on the sample, the collected data were normalized to the intensity of the incident 
beam. A typical IXS spectrum is shown on a logarithmic scale in figure H The experimental 
data with their error bars are shown together with the best fits (full lines) using both the 
model function discussed before (bottom panel) and the same function with At=0 (top 
panel). For clarity, only the elastic and the Stokes part of the spectrum is shown. The 
insets provide a zoom of the elastic line region on a linear scale. It can be easily appreciated 
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T = 263 K 




Energy (meV) 

FIG. 2: IXS spectra of water at T—263 K and P~2 kbar. The figure reports the data (dots) normahzed 
to the total integrated intensity on a logarithmic scale emphasizing the inelastic part of the spectrum. The 
solid line represents the total fit result as explained in the text, while the other lines visualize the two 
contributions to the fitting model: the longitudinal and quasi-elastic one represented by a viscoelastic model 
(dashed line) and the transverse-like secondary peak accounted for by a DHO model function (dash-dotted 
line) 

that a fit using a viscoelastic model only (top panel), does not properly describe the IXS 
data, leading to a poor fit of both the elastic line and the inelastic part of the spectrum (the 
reduced results to be ^ 2, i. e. more than 10cr^2 larger than its expected value). The 
inclusion of a second inelastic excitation significantly improves the quality of the fit (reduced 

~ 1, 0"x2 ~ 0.1), as can be seen in the lower panel. For all fitted spectra, the value of 
remained within one standard deviation from the expected value. 

Figure 121 shows as an example the Q evolution of the IXS data recorded at T=263 K and 
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T = 419K 




Energy (meV) 

FIG. 3: IXS spectra at r=419 K and P—95 bar. The figure reports the data (dots) normahzed to the total 
integrated intensity on a logarithmic scale emphasizing the inelastic part of the spectrum. The solid line 
represents the fit result as explained in the text, while the other lines visualize the two contributions to the 
fitting model: the longitudinal and quasi-elastic one represented by a viscoelastic model (dashed line) and 
the transverse-like secondary peak accounted for by a DHO model function (dash-dotted line) 

P=2000 bar. The two inelastic contributions which clearly appear at Q=10 nm^^ can also 
be seen throughout the whole Q-range. As can be noticed by the Q evolution of the dashed 
and dotted lines (longitudinal and transverse-like contributions, respectively) the intensity 
of the low-frequency (transverse-like) mode increases with increasing Q. 

The shape of the IXS spectra as well as their Q-dependence is markedly different at 
T=419 K and P=95 bar (see figure El). The elastic line appears broader, and consequently 
the low frequency, non dispersing mode is much less visible. As a matter of fact, it can 
only be identified in the spectra with Q < 13 nm^^, while for larger Q values the increasing 
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Energy (meV) 

FIG. 4: Temperature evolution of the IXS spectra for Q=10 nm~^ (left panel) and Q = 13 nm~^ (right 
panel). The figure reports the data (dots) normalized to the total integrated intensity on a logarithmic 
scale emphasizing the inelastic part of the spectrum. The solid line represents the fit result as explained 
in the text, while the other lines visualize the two contributions to the fitting model: the longitudinal and 
quasi-elastic one represented by a viscoelastic model (dashed line) and the transverse-like secondary peak 
accounted for by a DHO model function (dash-dotted line). 

width of the elastic hne completely governs the spectral shape. The weak feature at about 
45 meV in the spectra at Q=4 nm~^ is the contribution from the high-pressure cell windows 
and corresponds to the diamond longitudinal acoustic phonon. 

To show an example of the temperature evolution of the IXS spectra, and more specifi- 
cally, of the behavior of the second weakly dispersing feature, the spectra recorded at Q=10 
(left panel) and 13 nm~^ (right panel) are shown in figured at the indicated temperatures. 
At the lowest temperature, the experimental data are shown together with the best fit to 
the model function and the two individual components (viscoelastic model and DHO). For 
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FIG. 5: Q-dependence of the (scaled) structural relaxation time r* = x ri{T — 263K)/ri{T) (see text). 



The dashed line reports the fit results of 
valid up to Q = 7 nra~^ . 



1^ 



for the Q dependence of the relaxation time at T = 263-Rr, 



the higher temperatures, only the viscoelastic contribution and the total fit are shown for 
the sake of clarity and in order to emphasize the decreasing importance of the second exci- 
tation. We notice a clear trend for both momentum transfers: with increasing temperature 
the contribution of the second excitation becomes smaller and smaller to the extent that it 
completely disappears at the highest temperature. 



III. DATA ANALYSIS AND DISCUSSION 



A. The longitudinal dynamics 

In the present section we present and discuss the results obtained for the longitudinal 
dynamics which are retrieved from the viscoelastic analysis as described in detail in section 
II B. 

In figure El we report the Q-dependence of the structural relaxation time Tq, at different 
temperatures T, scaled by the shear viscosity ratio r(T) = rj{T = 263K)/ri{T): i. e. r* 
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FIG. 6: Temperature evolution of the structural relaxation time (full symbols) at different Q- values between 
7 and 16 rim~^. The shear viscosity 'ri{T) is reported for comparison (open dots), as well as the Q — 
extrapolation of the relaxation time 15] (solid line). 

= r{T)Ta{T). This scaling has been performed in order to verify whether the structural 
relaxation time is proportional to the longitudinal viscosity rj, a relationship which generally 
holds in the continuum limit Q — > 0. It is evident from the inspection of figure El that the Q 
evolution of at different T is indeed the same, especially in the Q-range below 20 nm~^. 
As implied by the validity of the scaling, the T and the Q dependence of the relaxation 
time can be factorized: Ta{Q, T) = ri(T)f{Q), being f{Q) the function that describes the Q- 
evolution of the relaxation time. The dashed line at low Q in Fig. reports the dependence 
of the relaxation time for T=263 K, derived in Ref. it turns out to be consistent with 
the present results. Figure El shows the temperature dependence of Tq, for different Q- values 
between 7 and 16 nm^^. The evolution of the viscosity (open circles) and the Q ^ limit 
of r„(0,T) (solid line) as reported in Ref. are also shown. The good overall agreemeirt 
confirms the reliability of the structural relaxation time determination by the viscoelastic 
model over the large Q range and in the different thermodynamic conditions investigated 
here. 

The temperature dependence of the zero frequency sound velocity (cq = uq/Q), the 
infinite frequency sound velocity (coo = \/C^M-u^/Q) and the apparent sound velocity 
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(ci = Ql/Q) are reported in figure [7| for T=263 K (left panel) and T=419 K (right panel). 
Here, fi^ is the maximum of the function iu'^s^^\Q, cu) calculated using the best fit parame- 
ters. The bottom panels of Fig.[7|show the Q-dependence of the product flLTa, the parameter 
that indicates whether the dynamics is viscous-like {^liTa < 1) oi elastic-like (^lTo > !)• 
For T=263 K, Qi'^a is always larger than, or close to, one, the system has a solid-like re- 
sponse, and the sound velocity cl is close to Coo, its infinite frequency value, over the whole 
explored Q-range. In contrast, at T=419 K, QlTo is larger than one only over a limited 
Q-range, from 4 to 10 nm~^, and then decreases rapidly below one for increasing Q- values. 
Consequently, we observe a transition from the infinite- to the zero-frequency sound regime: 
it takes place for Q values between 10 and 15 nm~^. This is the first experimental observa- 
tion of the transition Cqo — > Cq that takes place for Q-values around the Q's where the S{Q) 
shows its first maximum, as a consequence of the De Gennes narrowing. The Q dependence 
of the excitation frequency, in fact, shows a decrease with a minimum just in correspondence 
with the first maximum of the S{Q) as a memory of a "pseudo-Second Brillouin zone". At 
the same time the relaxation time decreases towards high Q with respect to the Q = value 
with a slight increase at around 27 nm~^ again due to De Gennes narrowing (see figure 
El). Such kind of Cqo — ^ Cq transition, taking place at Q values around the first maximum 
of the S{Q), has been already observed in MD simulations of a Lennard- Jones model glass 



31[. Furthermore, we note that our derived values for cq is in excellent agreement with 
independent determinations, both in the limit of Q ^ and at hig h Q- values. The values of 
Co in the low Q-limit were obtained from thermodynamic data j2l|, and are indicated by the 
arrows in Fig. [3 The zero frequency sound velocity in the high-Q region (dash-dotted line) 
has been calculated utilizing the expression for Cq{Q) within the framework of generalized 
hydrodynamics: 



Here, we utilized the S{Q) values determined by neutron diffraction measurements 
at the same thermodynamical conditions (the Sqo partial structure factor has been used). 
The agreement of co{Q), calculated using the previous relation, and cq{Q) derived from the 
fits convincingly demonstrates the solidity of our data analysis procedure. 

As expected for a viscoelastic behavior of the excitation, the same Cq ^ Cqo transition 
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can be observed at fixed Q due to the effect of the temperature on Tq,. Figure |H1 shows cl, 
Co and function of temperature for Q=1Q nm ^. At low T the relaxation time is 

long, VLiTa >> 1, and the apparent sound velocity is close to Coo- When, on increasing the 
temperature, the structural relaxation enters in the excitation timescale {VL^Ta ~ 1), the 
sound velocity cl recovers its liquid-like value Cq. 



B. The second excitation 



As already discussed in some detail above, the IXS spectra, especially at low temperature, 
can only be properly described if one takes into account a second inelastic excitation. This 
second excitation is described by three parameters (intensity, peak position and peak width). 
Figure ini shows the dispersion (Q-dependence) of VLt-, together with that of VLl, at T=263 K 
and P=2000 bar, where the second excitation is best appreciated. The upper panel shows 
the result of the present work. The longitudinal frequency VLl (squares) and the "transverse" 
frequency VLt (circles) (both derived from the fit procedure) show the expected behavior: a 
linear dispersion followed by a bend down and a minimum in the Q region where S{Q) has 
its maximum for the longitudinal excitation, and a weak Q dependence for the transverse 
mode. Full symbols in the lower panel correspond to the results of a previous MD simulation 
[isj l for the position of the Longitudinal current peaks. These were carried out considering 
4000 D2O molecules enclosed in a cubic box with periodic boundary conditions, and utilizing 
the SPC/E model 33|. The simulations were performed at a density of 1 g/cw? and T ^ 
250 K. The qualitative agreement between experiment and simulation is remarkable. We 
note slight differences in the absolute values of the excitation energies: while the maximum 
of the longitudinal dispersion is slightly lower, the energies for the second excitation are 
slightly higher in the case of the IXS results. Moreover a small displacement at higher Q 
in the position of the first minimum in the IXS data is justified as these are taken at high 
pressure while the simulations are performed at ambient conditions: as the pressure doesn't 
affect very much the value of Cqo the agreement with the simulation data in the linear region 
is extremely good. The open symbols, also reported in the lower panel, indicate for some 
Q values the position of the two peaks found in the Transverse Current spectra, a quantity 
that is not experimentally observable, but that can be determined by MD simulations. 
The position of these peaks clearly coincides with the correspondent peak measured in the 
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FIG. 7: Sound velocities and the product ^IlTu -as determined by the fit to the IXS spectra- as a function 
of Q, and at two different temperatures: T=263 K (left panels) and T=419 K (right panels). The top panels 
show the Q-evolution of the sound velocities cq (open dots), Coo (open triangles), and (full squares), 
whereas the bottom panels show the product ^lTu, that, when equal to unity, marks the transition from 
the infinite frequency sound velocity c^c to the zero- frequency sound velocity cq. The solid line is a guide to 
the eye. The hydrodynamic (Q=0) value of cq is indicated by arrows; the Q-evolution of co{Q), as derived 
using Eq. Eland the S{Q) determined by neutron scattering, is indicated by the dash-dotted line in the Q 
region where neutron data are available. The dotted horizontal line in the lower panels marks the condition 
^LTa = 1- 

Longitudinal Current, supporting the presence of a mixing phenomenon We further 

note that the energy range over which the weakly dispersing feature is observed corresponds 
to the TA / TO branch in hexagonal ice as determined by INS Q and IXS Q. 

To further characterize the second excitation we now study its integrated intensity. Specif- 
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FIG. 8: The temperature dependence of the sound velocities and of the product flLTa for Q=16 nm^^. 
The data show the same phenomenology as those reported in fig. |7| Here the product JIlTq is changed 
by varying the relaxation time scale: as in the previous case we start to observe the transition when we 
approach QlTu ~ 1- The solid line is a guide to the eye. The stars in the upper panel indicate the Co{Q) 
values as derived from Eq.|Hlfor the thermodynamical points available from the neutron diffraction database 

y 

ically, in figure ITUl we report the "transverse" intensity Aj- normalized to the total integrated 
intensity At + Al as a function of Q for the selected thermodynamic point T=263 K and 
P=2 kbar. In the lower panel of Fig. we report the evolution of Qt'^o- We observe a 
general increasing trend of the parameter At /{At + Al) with increasing Q. It is worth to 
observe that the appearance of a transverse-like contribution with an intensity increasing 
with Q has also been registered in experiments and simulations on glassy systems as glycerol 



3511 and silica 



27 



36 



37|. In silica the mixing phenomenon turns out to be particularly 
enhanced both in MD simulation and experiments indicating that, as in water, the local 
tetrahedral structure favors the coupling of the L and T dynamics. The intensity increase 
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FIG. 9: Panel (a): Longitudinal current peaks positions from the fits of IXS data of liquid water at T=263 



K and P—2 kbar. Panel (b): MD simulations results [18| (see text). The dotted lines represent excitations 
with a linear dispersion in Q, corresponding to a sound speed of 3300 and 1900 m/s, respectively. 

with Q is expected on the basis of a simple one-excitation picture within the harmonic 
approximation for the dynamic structure factor, suggesting that, for a non- dispersive exci- 
tation. At/ {At + Al) ^ Q^. More importantly, we also observe that, superimposed to the 
growth of At / {At + At) with Q, there is also a "modulation" in phase with Qt^o- When, at 
Q ~ 15 nm^^, I^tTq, become less than unity, the intensity of the "transverse" peak decreases. 

The presence of a strong correlation between the intensity of the transverse peak and the 
^^T^a is confirmed by inspecting their temperature evolution reported in figure ^2 for two 
selected Q values: Q =10 and 13 nm~^. Here, when the time scale of the a-relaxation is 
comparable to the inverse of the frequency of the transverse excitation, the ratio At/ {At + 
At) decreases, and for the highest temperature becomes essentially zero. This behavior -very 
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FIG. 10: Q-dependence of the transverse to total integrated intensity ratio (upper panel) and the product 
^TTa (lower panel) at T=263K. The dashed horizontal line indicates the condition QtTu = 1- 



different from that of tlie longitudinal collective mode, which is affected by the relaxation 
only in the sound velocity value- is just the one expected for a transverse-like excitation: 
recovering a liquid-like regime the system is no longer able to give an elastic response to 
a shear stress and to sustain propagating transverse waves. In this regime the transverse 
dynamics assumes a purely relaxational behavior, corresponding to a peak at = in the 
current spectrum . For this reason, still lacking a formal model for the coupling, the fits 
at higher temperatures have been performed just fixing the DHO parameters at the lowest 
T values leaving the intensity free. This procedure allows us to provide the information 
on the inelastic contribution of the dispersion-less mode leaving the relaxational behavior 
accounted for by the viscoelastic model. This procedure, as can be seen by bare eye in the 
spectra at high T, doesn't affect at all the statistical significance of the fits also when the 
transverse-like contribution, which was mandatory at low T, can be completely neglected. 



22 



Q=10nm"' Q=13nm"' 



a 



5- 

H 3- 



2- 

1- 

0- 
250 



300 



350 



400 250 



300 



350 



400 



T(K) 



FIG. 11: Top panels: temperature evolution of the transverse to total integrated intensity ratio At /{At + 
Al) at momentum transfers of 10 nm~^ (left panel) and 13 nm~^ (right panel). The solid lines are guides 
to the eye. Bottom panels: temperature evolution of OtTo- The dotted horizontal line marks the condition 
f^TTa = 1- 

IV. CONCLUSIONS 



In this paper the viscoelastic analysis of the IXS spectra of water, already successfully 
performed in the low Q range, has been extended to much larger Q values and in a wide 
range of thermodynamical conditions. 

At the high Q values investigated here the presence of a second excitation in the spectra 
cannot anymore be neglected and it has been successfully taken into account in the analysis 
revealing a behavior in agreement with the viscoelastic expectation for a transverse-like 
excitation. This findings confirm our assignment of a transverse nature to the dispersion- 
less mode of water observed at an energy of fa 6 meV. 
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Furthermore the large variations of the structural relaxation time (almost a decade) in 
the investigated thermodynamic range has allowed us to perform an excellent test of the 
viscoelastic picture. Summarizing, the two main results of the present work are: 

(1) The longitudinal-like dynamics can be properly accounted for by a viscoelastic model 
also in the high Q range: indeed, besides the now well known transition from cq to Coo 
that takes place at small Q (Q ~ 4 nm~^ at ambient conditions), for the first time the 
transition from Cqo back to Cq, taking place at high Q values, could be observed. This 
back-transition is a consequence of the De-Gennes narrowing on the product VtiTa 
that becomes smaller than one at Q values around the Q's of the main peak in the 
S{Q). 

(2) The second dispersion-less mode behaves as a transverse-like excitation disappearing 
from the spectra as soon as the structural relaxation process reaches the excitation 
time scale (fir'^a = 1) 

Finally, a few words must be devoted to the comparison between the present -viscoelastic- 
explanation of the high frequency dynamics of water and the alternative -solid based- de- 
scription as discussed in the introduction. To this purpose we compare in figure IT^ the tem- 
perature behavior of the longitudinal and transverse dynamics for a fixed Q value ((5=10 
nm~^). In panel (a) we report the velocities plot as in figure |H1 (here for a different Q value) 
while in panel (b) we report the (normalized) intensity of the dispersion-less mode at the 
same Q value. The comparison of the two uppermost panels in figure IT^ clear Iv demonstrates 
that while the apparent sound velocity (c^, full squares in panel (a)) is still far from from cq 
(open dots) -i. e. the system is in the elastic limit over the whole T range-, the intensity of 
the second excitation becomes negligible at high T. Thus we have that, at high T the system 
is in the elastic regime, but the second excitation can no longer be appreciated in the spectra. 
This evidence cannot be framed in the interaction model proposed in Ref. that traces 

the "fast sound" phenomenon back to the presence and the interaction between two modes. 
This observation, however, has a clear explanation in the viscoelastic picture as illustrated 
in panel (c). While the structural relaxation process has already reached the transverse-like 
mode time scales, with VLtToi < 1 for T > 370 K thus entering the viscous regime, it is 
still not affecting the higher energy longitudinal-like excitation which still exhibits an elastic 
response (fiiTc > 1 up to the highest explored temperature). 
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FIG. 12: Comparison of the temperature evolution of the L-like and T-like dynamics for a fixed Q value 
(Q= 10 nm^^). (a) Evolution of the three velocities as in figure |H1 the dashed line corresponds to the 
thermodynamic prediction for cq. (b) Temperature evolution of the transverse to total integrated intensity 
ratio At /{At + Al). (c) Temperature evolution of JItTq and QlTu] the dotted horizontal line marks the 
condition flTa ~ 1 clearly showing that the two dynamics are in different regimes at high T. The solid lines 
are guides to the eye. 

The data reported in Fig. together with the outcome hsted in items (1) and (2) 
before, are the main results of this paper. All these results favor the viscoelastic description 
of the high frequency dynamics of water. As a final observation, it must be noted that 
the two models discussed before are not so different as it appears at a first glance: in both 
cases it is the interaction of the sound waves with another "mode" which causes both the 
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positive dispersion of the apparent sound velocity and the appearance of a second peak in 
the dynamic structure factor. In the viscoelastic case the second mode is purely relaxing, 
in the solid-like approaches the second modes have a propagating component. However, 
major differences arise in interpreting the physical meaning of the parameters entering in 
the models: i) in the solid-like model no information of the structural relaxation time is 
contained, while is directly measurable in the viscoelastic model, it is in agreement with 



independent determinations 



15l |. and its values allow to calculate the viscosity, again 
in agreement with independent measurements ii) A further evidence, again in the 

same direction, comes from the temperature dependence of the energy of the second mode. 
In the solid-like framework, at increasing temperature the data indicate a decrease of the 
interaction parameter, corresponding to a decrease of the energy position of the dispersion- 
less mode. This decrease of the energy of the dispersion-less mode is not observed; for 
example, the depolarized light scattering measurements js^ do not show any noticeable shift 
in the position of the peak that corresponds to the density of states of the second excitation; 
Hi) Finally, let us observe that the viscoelastic-based explanation of the dynamics of the 
de,.,ty fluctuations . actually considered the proper approach in many l„u,ds Q and rt 



has been substantiated by rigorous theories as the Mode Coupling 



heory j40| (a theory 



). In other words. 
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that, via MD, has been proved to work properly also in liquid water 
as the dynamic of liquid water shows the same phenomenology of that of many other liquids, 
why, in water, should one look for a different explanation? 
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